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Abstract 

For  laminar,  steady  flow  in  ducts,  the  current 
definition  of  hydraulic  diameter,  Dk,  does  not  accurately 
depict  the  non-uniform  wall  shear  stress  distribution  around 
the  perimeter  of  non-circular  duct  shapes.  In  this 
investigation,  a  new  hydraulic  diameter,  D„.i,  was 
empirically  determined.  It  correlated  friction  factor  data 
for  many  non-circular  shapes  to  within  approximately  2.4  % 
of  the  circular  duct  value. 

An  experiment,  using  the  AFIT  Oil  Flow  Rig  Set-up,  was 
run  to  determine  the  effect  on  transition  Reynolds  number. 
Re.,,  and  hydrodynamic  entrance  length,  L*,  of  replacing  Dh 
with  Dh.i.  Transition  Reynolds  number  and  L*  were 
determined,  based  on  Dh  and  Dk,»,  for  a  circular,  square, 
and  concentric  annular  duct. 

Transition  Reynolds  numbers,  based  on  Dk,  for  the 
square  and  concentric  annular  ducts  were  approximately 
12.5  %  lower  than  the  circular  duct  Re,,.  The  Re.,,  based 
on  Dh, i ,  did  not  correlate  well  for  the  concentric  annulus, 
but  did  correlate  for  the  square  duct. 

Hydrodynamic  entrance  lengths,  based  on  Dk  and  Dh,i, 
were  experimentally  determined  for  the  circular  and  square 
duct  only.  The  square  duct  L*,  and  the  analytic  concentric 
annular  duct  L*,  based  on  Dk.i,  did  not  correlate  to  the 


to  the  circular  duct  value. 

Although  the  Re*,  and  L*  for  the  square  and  concentric 
annular  ducts  did  not  correlate  well  when  based  on  Dk.i,  th< 
data  obtained  is  still  useful  to  the  engineering  community 
since  it  provides  an  addition  data  base  for  experimental  L 
and  provides  experimental  data  on  Re*,  for  square  and 
concentric  annular  ducts. 


A  REDEFINED  HYDRAULIC  DIAMETER  FOR  LAMINAR  FLOW 


I .  Introduction 


Background 

For  flow  in  circular  ducts,  where  there  is  a  uniform 
wall  shear  stress  around  the  perimeter  and  uniform  distance 
from  the  center  to  the  surface,  duct  radius,  or  diameter, 
is  a  logical  choice  for  the  characteristic  dimension  in 
flow  correlations.  For  non-circular  ducts  a  characteristic 
dimension,  the  hydraulic  diameter  (D„),  has  been  defined  as 
4 A/P  where  A  is  the  duct  cross-sectional  area  and  P  is 
the  duct  perimeter.  The  hydraulic  diameter  is 
unambiguously  defined  for  any  cross-sectional  shape,  and 
for  consistency  is  used  in  both  turbulent  and  laminar  duct 
flows.  Experimentation  has  verified  that,  for  turbulent 
flow,  where  the  wall  shearing  stress  is  fairly  uniform 
around  the  perimeter,  the  hydraulic  diameter  is  the 
appropriate  characteristic  dimension  for  many  non-circular 
duct  shapes  (1).  For  fully  developed,  laminar  flow, 
however,  the  wall  shearing  stress  distribution  along  the 
perimeter  of  non-circular  ducts  is  not  uniform  and  the 
hydraulic  diameter  is,  perhaps,  not  the  best  characteristic 
dimension.  For  example,  for  a  circular  duct,  C»Re  =  16.0, 
while  for  a  triangular  duct,  C,Re  =  13.333,  where  Re  is 
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defined  as  VD*/v.  At  present  there  is  no  characteristic 
dimension  that  correlates  laminar  flow  data  in  non-circular 
ducts . 

Some  current  uses  of  hydraulic  diameter  as  a 
characteristic  dimension  are  in  transition  Reynolds  number 
and  hydrodynamic  entrance  length  calculations.  Reynolds 
number,  which  is  a  ratio  of  inertia  to  viscous  forces,  is 
very  much  associated  with  the  stability  of  laminar  flows. 
Viscous  forces  tend  to  restore  laminar  flow  after  a 
disturbance,  while  inertia  forces  tend  to  amplify 
disturbances, therefore,  transition  from  laminar  to 
turbulent  flow  might  be  characterized  by  a  transition 
Reynolds  number,  Re.,.  For  circular  duct  flow,  the  Re.,, 
based  on  diameter  and  duct  averaged  velocity,  is 
approximately  2000  (2).  This  author  has  not  found  any 
Re.,  data  for  non-circular  ducts,  and  since  the  wall 
shearing  stress  is  not  uniformly  distributed  around  the 
perimeter  of  non-circular  ducts,  one  would  speculate  that 
the  Re*,,  based  on  hydraulic  diameter,  for  these  ducts 
would  be  different  from  the  circular  duct  value  of  2000. 

Another  use  of  hydraulic  diameter  has  been  in 
hydrodynamic  entrance  length  (X)  calculations.  One  of  the 
most  common  definitions  of  X  is  the  duct  length  required  to 
achieve  a  duct  centerline  velocity  99%  of  the  fully 
developed  value  (3).  These  results  are  often  given  in  the 
non-dimensional  form,  L*  =  (X/D».)/Re.  Several  different 
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analytical  solutions  to  the  L*  problem  exist  for  circular 
as  well  as  non-circular  duct  shapes  (14).  Based  on 
hydraulic  diameter,  L*  varies  significantly  for  each  duct 
shape.  Little,  if  any,  experimental  L*  data  has  been 
reported. 

Objectives 

This  investigation  had  four  objectives: 

1.  Empirically  determine  a  new  hydraulic  diameter, 
"Dh.i " ,  that  better  represents  the 

mean  center  to  surface  distance  for  laminar 
flow  in  non-circular  ducts,  and  which 
correlates  CfRe  data. 

2.  Experimentally  determine  the  Re*,  ,  based  on 
Du,* ,  for  a  circular  duct,  square  duct,  and 
concentric  annulus  with  an  aspect  ratio  that 
simulates  flow  between  infinite  parallel 
plates . 

3.  Experimentally  determine  L*  for  the  above 
shapes,  based  on  D„.i. 

4.  Examine  the  effect  on  Re*,  and  L*  of 
replacing  the  traditional  D*  with  Dh.i 

The  first  objective  will  be  accomplished  using 
friction  factors  (CfRe),  for  fully  developed  duct  flow,  as 
the  data  base  for  the  Dh.i  calculations,  with  an  attempt  to 


converge  the  data  to  the  circular  duct  value. 

The  remaining  objectives  will  be  accomplished  using 


the  Air  Force  Institute  of  Technology's  instrumented  oil 
flow  rig  set-up.  Pressure  and  mass  flow  data  will  be  taken 
for  the  different  duct  shapes  to  determine  Re.r  and  L*. 


C? 


I I .  Duct  Theory 


Circular  Ducts 

The  circular  duct  is  the  most  common  duct  shape  in  use 
today.  For  constant-property,  steady,  fully-developed, 
laminar  flow  through  circular  ducts,  as  in  Figure  1,  the 
equation  of  motion  can  be  written  as: 

(M/r ) [ <d/3r ) (r(*u/dr ) ) ]  =  dP/dx  (1) 

Since  the  pressure  is  independent  of  r.  Equation  (1)  can 
be  integrated  directly  with  respect  to  r  to  obtain  the 
velocity  profile.  Applying  boundary  conditions: 


r  =  0 


axis  of  the  tube 


Su/ar  =  0 


velocity  gradient  is  zero 


Z  -  Zo 


duct  wall 


u  *  0 


no  slip  at  the  wall 


yields 


u  =  (ZoV4n)  (-dP/dx)  (l-(r*/ro2))  (2) 

It  proves  more  convienient  to  express  this  velocity  in 
terms  of  a  mean  duct  velocity,  V,  rather  than  the  pressure 


gradient.  For  incompressible  flow,  the  mean  velocity  13 


defined  as: 


resulting  in. 


thus. 


-«// 


udA) /A 


V  =  (ro*/8M)  (-dP/dx) 


u  =  2V(  1-  ( z/io)3) 


Entrance 


Entrance 


ruiiy 


Region 


Developed 


Figure  1.  Flow  in  a  Circular  Duct 


The  gradient  of  this  velocity  profile  is  used  to  evaluate 


the  shear  stress  at  the  wall. 


r*  ■  m  (du/dr  ).-»• 


Using  Equation  (4) 


%  »  Ml2V(-2r/r.*)  1^*  -4VM/r. 


Since  V  is  constant  for  steady,  fully-developed  flow,  u  is 


constant.  The  friction  factor  can  now  be  evaluated.  The 


Fanning  friction  factor  (C»)  is  defined  as  the  ratio  of 


wall  shearing  stress  to  the  dynamic  head  of  the  flow. 


C,  =  %/(PV*/2) 


thus , 


C,  =  ( 4V<i/r«)/(PV*/2 )  =  8M/PVro 


'  '  ■'  W'nV-.V-.V 


For  a  circular  duct.  Re  is  defined  as 

Re  =  PVD/fx 

where  0  =  2r<>.  Rewriting  Equation  (8) 

C,  =  16/Re  (9) 

Data  from  Stanton  and  Pannell  (6:199)  and  Senecal  and 
Rothfus  (7:533)  verifies  this  equation  up  to  a  Re  of  2000. 


Won-Circular  Ducts. 


Although  a  majority  of  ducts  in  use  today  are  of 
circular  cross-section,  many  non-circular  duct  applications 
exist  and  require  further  flow  analysis.  For  constant 
property,  fully  developed,  laminar  flow  through  a 
non-circular  duct,  the  equation  of  motion  is: 

3au/3y*  ♦  3au/aza  =  (1/fi)  OP/3x) 

The  friction  forces  are  seen  to  be  functions  of  the 
cross-sectional  dimensions  in  both  y  and  z  directions. 
Traditionally,  a  single  length  dimension,  the  hydraulic 
diameter  (D»),  has  been  used  to  characterize  the 
cross-section .  Hydraulic  diameter  is  related  to  the 
hydraulic  radius,  R»,,  which  is  the  ratio  of  cross-sectional 
area.  A,  to  wetted  perimeter.  Per,  by: 

D*  =  4R*  =  4  A/Per 

Hydraulic  diameter  for  a  non-circular  duct  is  intended  to 
represent  the  diameter  of  an  equivalent  circular  duct  which 
has  the  same  C> Re  as  the  non-circular  duct.  Table  I 


illustrates  D*  for  several  non-circular  shapes. 


TABLE  I. 


Five  Ducts  With  the  Same  Hydraulic  Diameter 


Oon  vctiOn 


Circular 


Square 


Infinite  paiallt: 
plates 


.Dh=2s  | 


5:77) 


Turbulent  flow  through  non-circular  ducts,  using  D» 


as  the  characteristic  dimension,  was  investigated  by 
Shiller  (8)  and  Nikuradse  (9),  who  verified  that  Ok 
provided  satisfactory  results.  For  laminar  flows,  however, 
the  Shiller  and  Nikuradse  data,  as  well  as  that  of  Koch  and 
Feind  (10),  showed  limitations  to  the  use  of  hydraulic 
diameter.  This  data  is  shown  in  Figure  2  from  Schlicting 
(2:576).  Here  A  is  the  Blasius,  or  Darcy,  friction  factor 
where, 

A  *  64/Re  =  4Cf 

for  circular  ducts.  In  Figure  2,  A  has  been  multiplied  by 
a  constant,  for  each  duct  shape,  in  order  to  fit  all  the 
data  onto  one  graph.  For  laminar  flow,  the  hashed  lines 
represent  a  circular  duct.  Note  the  difference  in  A 
between  circular  and  non-circular  ducts  based  on  D». 
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Figure  2.  Friction  Factor  for  Several  Ducts 


Exact  solutions  to  the  equations  o£  motion,  for  many 
non-circular  ducts,  based  on  D*.,  yield, 

C»Re  =  C  =  constant 

for  each  duct  shape,  instead  of  the  desired  circular  duct 
solution,  C,Re  =  16.  Several  values  of  C*Re  are  shown  in 
Table  II  (11). 

TABLE  II. 

CrRe  for  Non-Circular  Ducts 

DUCT  SHAPE  C  =  C>Re 


SINE 

12.630 

TRIANGLE 

13.333 

SQUARE 

14.227 

HEXAGON 

15.054 

RECTANGLE  2:1 

15.548 

SEMI -CIRCLE 

15.600 

RECTANGLE  4:1 

18.233 

RECTANGLE  8:1 

20.585 

INFINITE  PARALLEL 

PLATES 

24.000 

These  values  of  C,Re  represent  exact  solutions,  based  on 
D*,  most  having  been  verified  by  experimental  data 
(8,9,10) . 


By  examining  the  definition  of  hydraulic  radius,  an 


insight  as  to  why  these  laminar  flow  values  do  not 
correlate  with  the  equivalent  circular  duct  data  can  be 
gained.  The  reciprocal  of  hydraulic  radius  is  the  wetted 
perimeter  per  unit  of  flow  cross-section  and  describes  the 
amount  of  wall  area  in  contact  with  the  moving  fluid.  The 
hydraulic  radius  is  a  good  index  if  the  resistance  to 
momentum  transfer  is  predominantly  dependent  on  the  amount 
of  wall  surface  area,  as  it  is  in  turbulent  flow.  In 
laminar  flows  this  resistance  to  momentum  transfer  takes 
place  throughout  the  cross-sectional  flow  area  and  is  not 
primarily  associated  with  the  wall  surfaces.  Thus  reliable 
conversions  from  circular  to  non-circular  shapes  cannot  be 
expected . 

Another  explanation  for  poor  correlation  with  laminar 
flows  is  seen  by  examining  the  surface  area  per  unit  of 
flow  volume,  As/Vol,  for  individual  ducts.  The  friction 
factor  is  directly  related  to  the  amount  of  wall  shear 
stress  and  this  stress  is  related  to  the  amount  of  surface 
area  it  acts  upon.  For  ducts  with  As/Vol  different  than 
that  of  a  circular  duct,  the  shearing  stress  is  acting  on  a 
different  amount  of  surface  area,  per  unit  volume,  and  a 
different  Cf  can  be  expected.  For  example,  a  rectangular 
duct  with  an  aspect  ratio  of  8:1,  has  more  surface  area, 
per  unit  volume,  for  the  shearing  stress  to  act  upon  than 
does  the  circular  duct  with  the  same  Dh,  and  hence  a  larger 


C*  is  expected. 


Flow  in  the  corners  of  non-circular  ducts  also 
influences  C. .  Consider  a  square  duct  and  a  circular  duct, 
each  with  the  same  Dh.  Each  duct  has  approximately  the 
same  As/Vol  ratio.  In  the  corners  of  the  square  duct,  the 
velocity  is  low  compared  to  V,  and  with  this  low  velocity 
is  low  shearing  stress.  This  low  velocity,  in  essence, 
reduces  the  effective  flow  area.  Low  effective  flow  area 
produces  a  lower  C<  for  the  square  duct  compared  to  the 
circular  duct. 

Since  laminar  flow  C.Re  values  do  not  correlate  well 
for  non-circular  ducts,  it  becomes  desirable  to  have  exact 
solutions  for  non-circular  shapes,  based  on  a  new  Dh.i, 
where  the  value  of  C.Re  for  these  ducts  equals  that  of  the 
circular  duct.  Individual  C.Re,  as  in  Table  II,  would  no 
longer  be  necessary. 

Exact  C.Re  values  vary  non-linearly  for  rectangular 
ducts,  from  the  square  duct  value  of  14.227,  to  the 
infinite  parallel  plate  duct  value  of  24.0,  as  illustrated 
in  Figure  3  (1).  Since  the  corners  of  rectangular  ducts 
have  the  same  influence  on  C.,  from  duct  to  duct,  the  C.Re 
variance  can  be  directly  attributed  to  the  difference  in 
As/Vol.  Once  again,  the  values  in  Figure  3  are  based  on 
Dh .  As  noted  in  Figure  3,  a  rectangular  duct  with  an 
aspect  ratio  of  approximately  2.3  has  a  C.Re  equa]  to  that 
of  a  circular  duct,  and  Dh  becomes  an  appropriate 


Currently,  each  non-circular  duct  has  a  unique  exact 
solution  for  C,Re,  based  on  D*.  This  author  seeks  a  new 
parameter,  Dn.i,  that  would  produce  a  non-circular  duct 
value  of  CtRe  equal  to  the  circular  duct  value  of  16.  The 
first  step  is  to  examine  the  C.Re  data  for  several 
non-circular  duct  shapes,  where 

Re  =PV  Dh/M 

Table  II  and  Figure  3  summarize  C* Re,  based  on  D*,  for  many 
of  the  most  common  non-circular  shapes.  An  empirical 
equation  for  Di.,1  will  be  derived  that  will  correlate  C.Re 
data  to  the  circular  duct.  To  do  this,  first  the  data  from 
square,  triangular,  hexagonal,  and  semi-circular  ducts  are 
examined,  since  these  shapes  are  similiar,  in  that  they  are 
fairly  compact.  Next  the  rectangular  data  from  Figure  3  is 
approximated,  to  within  0.05%,  by  the  following  equation: 
CfRe  =  2  4 [  1  -  (  1.3553/a)  *  (  1.9467/a*)  -  (1.7012/a») 

+  (.9564/a4)  -  (.2537/a8)]  (3:199) 

The  rectangular  duct  will  be  examined  first  for  a  from  1.0 
to  10.,  in  increments  of  0.25.  Most  manufactured 
rectangular  ducts  are  in  this  range,  so  small  increments 
are  taken  to  fully  evaluate  these  ducts.  Aspect  ratios  of 
25,  50,  100,  and  »  were  examined  next.  These  ducts 


represent  flow  that  is  beginning  to  simulate  infinite 
parallel  plate  flow. 

The  hydraulic  diameter  describes  the  diameter  of  an 
equivalent  circular  duct  which  has  the  same  mean  shear 
stress  as  the  non-circular  duct.  Since  Dh  under-estimates 
the  correct  diameter  in  some  cases,  such  as  square  and 
triangular  ducts,  and  over-estimates  it  in  others,  such  as 
infinite  parallel  plates  and  elongated  rectangular  ducts,  a 
new  Oh,  Dh.!,  is  invisioned.  Dh,i  is  determined  using  a 
"modelling  circle"  diameter.  Dm,  plus  a  correction  factor. 
This  "modelling  circle"  represents  the  area  of  "main  flow" 
in  the  duct.  It,  in  essence,  separates  the  corner  areas 
from  the  "main  area"  in  the  duct.  Table  III  illustrates 
the  modelling  circles  for  several  duct  shapes.  For  fairly 
compact  shapes,  such  as  square,  triangular,  and  hexagonal 
ducts,  the  modelling  circle  area  is  identically  the  Dh 
area,  and  D.  =  Dh.  For  elongated  ducts,  such  as  infinite 
parallel  plates  or  rectangular  ducts,  where  the  flow 
extends  out,  the  modelling  cirle  is  an  ellipse  where  the 
elliptic  aspect  ratio  is  the  same  as  that  of  the  duct.  For 
simplicity,  D.  is  the  smaller  length,  a,  for  elliptical 
modelling  circles.  Comparing  the  flow  area  of  circles  with 
diameter  Dh  and  their  associated  perimeter  to  those  of  the 
elliptical  modelling  circles,  shows  that  the  elliptical 
modelling  circles  better  represents  the  actual  flow  area, 
and  the  perimeter  upon  which  the  mean  shear  stress  may  act. 


TABLE  III 


Modelling  Circles  for  Several  Duct  Shapes 


D.  «  1.000  S 
Dk  =  1.000  S 


D.  =  .57735  S 
Dk  *  .57735  S 


D.  *  1.732  S 
Dk  *  1.732  S 


2b /2a  - 


Dk  *  .500  S 
Dk  *  .611  S 


2b /2a 


D.  =  1.000  S 
Dk  =  1.333  S 


1.000  S 
1.600  S 


1.000  S 
2.000  S 


16 


This  can  be  seen,  for  example,  for  the  rectangular  duct  of 
a  =  4:1,  by  comparing  the  sketches  in  Table  I  and  Table 
III.  In  determining  a  correction  factor  for  D.,  the  area 
between  the  duct  walls  and  the  modelling  circle,  and  its 
associated  perimeter,  is  utilized. 

Results 

For  laminar  flow  in  non-circular  ducts: 

Dh..  =  3 . 2  { AE/PE  )  (  Dh/D, )  +  D.  (10) 

where 

AE  =  area  of  the  duct  -  area  of  the  "modelling  circle' 
PE  =  perimeter  of  the  duct  +  perimeter  of  the 
"modelling  circle" 

Dh  =  hydraulic  diameter,  4A/P 

D.  =  smallest  diameter  of  elliptic  "modelling  circle" 
Figure  4  shows  an  example  of  AE  and  its  associated 
perimeter.  Equation  (10)  is  used  for  all  ducts  except 
rectangular  ones.  An  additional  term  is  needed  for 
rectangular  ducts  to  account  for  the  changing  Dh/D.  ratio 
and  the  aspect  ratio,  a.  For  rectangular  ducts: 

D„,.  =  3. 2  (AE/PE)  (Dh/D-)  [1  +  (  4  ( Dh  -D.)-9/a)]  +  D.  (11) 

For  example,  using  Figure  4  and  Equation  (11), 
area  of  the  duct  =  3S* 

area  of  the  "modelling  circle"  =  * ( . 5S ) ( 1 . 5S ) 

perimeter  of  the  duct  =  8S 

perimeter  of  "modelling  circle"  =  2  jt[  ( 1 . 52  +  .5*)/2)-9S 


Figure  4.  AE  for  a  Rectangular  Duct 

AE  =  . 6  44S2 
PE  =  15.02S 
aspect  ratio  =  3.0 
D„  =  1.5S 
D.  =  1.0S 

and, 

D*.,  =  1.4S 

Here,  for  an  ellipse: 

Area  -  irab 

and  the  perimeter  can  be  approximated  by. 

Perimeter  2*1  (a2  +  b2)/2)-8 
Note  that  equation  (11)  reduces  to  equation  (10)  for  the 
limiting  cases  of  the  square  duct  ( Dk  =  D.)  and  the 
infinite  parallel  plate  duct  (a  -<*>). 

Table  IV  and  V  present  Da.,  and  C.Re  data,  based  on 
Da. i  .  Friction  factors  for  all  duct  shapes  investigated, 
are  within  2.4%  of  the  circular  duct  value  of  16.0. 
Appendix  A  presents  a  breakdown  of  the  Da.i  calculations 


for  each  duct  shape. 


SHAPE 

CfRe  (Dh) 

Dh 

:ffte  (Dh 

SQUARE 

14.227 

1.000S 

1.000S 

1.098S 

15.62 

TRIANGLE 

13.333 

0.577S 

0.577S 

0.693S 

15.94 

HEXAGON 

15.054 

1.732S 

1.732S 

1.800S 

15.64 

SEMI -CIR 

15.600 

0.611S 

0.500S 

0.633S 

16.15 

IHF.PAR.PL.  24.000  2.000S  1.000S  1.345S  16.12 
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TABLE  V. 


T*rr 


T^rx^i 


! 


Redefined  Hydraulic  Dlaroeter^Tiij-1 
Rectangular  Bucia 


Aaassl 

Ratio 

CfRe  (D^) 

1 1  1  — 

• 

1.00 

14.230 

1.000S 

1.098S 

15.620 

— 

2.4 

1.25 

14.382 

1.111S 

1.239S 

16.039 

0.2 

1.50 

14.715 

1.200S 

1.294S 

15.863 

0.9 

1.75 

15.127 

1.273S 

1.328S 

15.786 

— 

1.3 

2.00 

15.557 

1.333S 

1.352S 

15.770 

1.4 

2.25 

15.978 

1. 385S 

1.368S 

15.783 

- 

1.4 

2.50 

16.377 

1.429S 

1.379S 

15.810 

- 

1.2 

2.75 

16.749 

1. 467S 

1.387S 

15.843 

- 

1.0 

3.00 

17.095 

1.500S 

1.393S 

15.877 

— 

0.8 

3.25 

17.414 

1.529S 

1.397S 

15.909 

— 

0.6 

3.50 

17.710 

1.556S 

1.400S 

15.940 

- 

0.4 

3.75 

17.982 

1.579S 

1.402S 

15.969 

- 

0.2 

4.00 

18.234 

1.600 a 

1.404S 

15.995 

• 

0.0 

4.25 

18.467 

1.619S 

1.404S 

16.019 

0.1 

4.50 

18.684 

1 .6365 

1.405S 

16.040 

0.3 

4.75 

18.885 

1.652S 

1.405S 

16.060 

0.4 

5.00 

19.072 

1.667S 

1.405S 

16.078 

0.5 

5.25 

19.246 

1.680S 

1.405S 

16.095 

0.6 

5.50 

19.409 

1.692S 

1.405S 

16.110 

0.7 

5.75 

19.561 

1.704S 

1. 404S 

16.124 

0.8 

6.00 

19.705 

1.714S 

1.404S 

16.136 

0.9 

6.25 

19.839 

1.724S 

1.403S 

16.148 

0.9 

6.50 

19.965 

1.733S 

1. 403S 

16.159 

1.0 

6.75 

20.084 

1.742S 

1. 402S 

16.170 

1.1 

7.00 

20.197 

1.750S 

1.402S 

16.180 

1.1 

7.25 

20.303 

1.758S 

1. 401S 

16.189 

1 . 2 

7.50 

20.404 

1.765S 

1.401S 

16.197 

1.2 

7.75 

20.499 

1 . 77 IS 

1. 400S 

16.206 

1 .  3 

8.00 

20.590 

1.778S 

1. 400S 

16.214 

1 . 3 

8.25 

20.676 

1.784S 

1.399S 

16.221 

1 . 4 

8.50 

20.758 

1.789S 

1.399S 

16.229 

1 . 4 

8.75 

20.836 

1. 795S 

1. 399S 

16.236 

1  .  5 

9.00 

20.910 

1.800S 

1.398S 

16.242 

1.5 

9.25 

20.981 

1. 805S 

1. 398S 

16.249 

1.6 

9.50 

21.049 

1.810S 

1. 397S 

16.255 

1.6 

9.75 

21.114 

1.814S 

1 . 397S 

16.262 

1.6 

10.00 

21.176 

1 . 818S 

1. 397S 

16.268 

1.7 

25.00 

22.770 

1 .923S 

1.370S 

16.218 

1 . 4 

50.00 

23.370 

1.961S 

1. 358S 

16 . 178 

1 . 1 

100.00 

23.680 

1.980S 

1. 351S 

16 .154 

1.0 

00 

24.000 

2 . 000S 

1. 343S 

16.120 

0 . 8 
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The  AFIT  Oil  Flow  Rig  Set-Up,  shown  in  Figure  5,  was 
used  to  determine  Re.,,  L*,  and  C<  for  a  circular,  square, 
and  concentric  annular  duct.  Oil  is  pumped  from  the 
reservoir,  up  to  the  duct  entrance,  down  the  duct,  and 
either  into  a  weighing  tank  or  back  into  the  reservoir, 
depending  on  the  selector  valve  setting.  A  mercury 
manometer  board  displays  individual  pressures  for  10  static 
pressure  stations  along  the  duct. 

SAE  10W-10  oil  is  the  working  fluid  for  this 
experiment.  To  determine  Re.,,  L*,  and  C,,  the  density  (  p) 
and  the  viscosity  (p)  of  the  oil  are  needed.  The  density 
of  the  oil  is  determined  by, 

p..i  =  Pmofspecific  gravity  of  oil) 

The  density  of  water  is  62.4  lbm/ft*,  and  the  specific 
gravity  of  oil  is  the  ratio  of  the  mass  of  a  given  volume 
of  oil  to  that  of  an  equal  volume  of  water.  Six 
measurements  were  taken  with  an  average  specific  gravity, 
for  this  oil,  of  .863  (see  Appendix  B.).  Thus,  the  density 
of  the  oil  is  53.85  lbm/ft*. 

The  viscosity  of  the  oil  can  be  determined  by  using 
the  Hagen-Poiseuille  law  (20:114),  or  by  use  of  a 
viscometer.  Lieutenant  J.  C.  Ghiglieri  (29),  determined 


Fi (lire  5 •  Oil  Flow  Rip  Set-up 


100  F 


10.18  centistokes 


210  F  -  2.53  centistokes 

Using  a  Viscosity-Temperature  chart,  like  that  shown  in 
Figure  6,  these  values  can  be  linearized  and  p  at  other 
temperatures  extrapolated.  Table  VI  shows  some 
extrapolated  values.  These  values  were  verified  from  the 
circular  duct  data  using  the  Hagen-Poiseuille  law, 

dP/dx  =  (128  n  Q  ) /  (  n  D4) 

where, 

Q  =  m/p  =  AV 

Circular  Duct 

Figure  7  shows  a  diagram  of  the  circular  duct  that  was 
used  for  this  experiment.  The  duct  is  21  feet  long,  with 
an  inside  diameter  of  0.8125  inches.  Ten  static  pressure 
taps  are  mounted  along  the  duct  with  the  first  tap  8  inches 
from  the  entrance,  and  the  remaining  9  taps  located  every  2 
feet  downstream.  The  cross-sectional  flow  area  is  0.5185 
square  inches. 

Square  Duct 

The  square  duct,  shown  in  Figure  8,  is  21  feet  long, 
with  an  inside  width  of  0.75  inches.  Pressure  taps  are 
mounted  identical  to  those  of  the  circular  duct.  The 
cross-sectional  flow  area  for  this  duct  is  0.5625  square 


inches . 


Tefflperoitf't.  ot9  f 


V lscos  1  ty-T emperature  Chart  Showing  Typical  Experimental  Curves 
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TABLE  VI . 

AFIT  SAE  10W10  Oil  Viscosity 


TemD. -F 

v  -  Centistokes 

U -  ( lbf-sec ) /f t 

70 

18.80 

3.39  X  10- 

75 

17 . 50 

3.15 

80 

15.75 

2 .84 

85 

13.92 

2.51 

90 

12.86 

2.32 

95 

12.07 

2.18 

100 

10.18 

1.83 

105 

9 . 50 

1.71 

110 

9.00 

1.62 

115 

8 . 20 

1 .48 

120 

7 . 50 

1 .35 

125 

6.90 

1 .24 

25 


» 


r 

* 

: 


Entrance 


Scale  in. 


A  -  .5185  in. ^ 

Figure  7.  Circular  Duct  Diagram 


Scale  in. 

A  =  .5625  in.2 


Figure  8.  Square  Duct  Diagram 


r,/r,  =  .66/. 7825  =  .843 5 

A  =  .533  in.2 


Figure  9.  Concentric  Annulus  Diagram 

Concentric  Annulus 

The  concentric  annulus,  shown  in  Figure  9,  consists  of 
two  ducts,  one  mounted  within  the  other,  with  r./r*  = 

.8435.  The  center  duct  is  blocked  at  both  ends  to  allow 
flow  between  the  two  ducts  only.  The  concentric  annulus 
was  designed  to  simulate  infinite  parallel  plates  flow, 
which  represents  the  upper  bound  of  C*Re .  Figure  10  shows 
that  infinite  parallel  plates  can  be  simulated  with 
concentric  annuli  having  r,/zi  greater  than  approximately 
0.7  (3:286).  This  duct  is  7.5  feet  long  with  a  flow  height 
of  .1225  inches.  This  results  in  a  cross-sectional  flow 
area  of  .555  square  inches.  Six  pressure  taps  are  mounted 
on  the  duct  with  the  first  three  taps  located  3.0,  5.3,  and 
8.0  inches,  respectively,  downstream  of  the  entrance.  the 
remaining  three  taps  are  located  every  two  feet  downstream 
from  tap  number  3.  Since  the  entrance  region  developes 


is  presented  on  the  mercury  manometer  board.  The  weighing 
tank  is  used  to  measure  the  mass  flow  of  oil  for  a  specific 
time  period.  By  dividing  the  total  mass  weighed  by  the 
elapsed  time,  a  mass  flow  rate  can  be  determined. 

m  =  lbm  oil/time 

For  friction  factor  measurements,  steady,  laminar  flow 
is  established  in  the  duct,  and  static  pressure  is  recorded 
for  each  of  the  10  stations,  from  the  manometer  board. 
dP/dx  is  determined  from  the  fully  developed  (linear) 
portion  of  the  pressure  curve  (see  Figure  11).  The  average 
velocity  (V),  is  determined  from  mass  flow  rate  data  where, 

V  =  m/p  A 

Friction  factors  are  then  determined  from, 

4Cf  =  (  -dP/dx )  (  2gD/P V* ) 

Transition  from  laminar  to  turbulent  flow  was 
determined  visually  at  the  duct  exit,  with  the  aid  of  a 
strobotac.  Once  laminar  flow  is  established  in  the  duct, 
the  flow  velocity  was  increased  until  turbulent  eddies 
become  noticable  at  the  duct  exit.  The  strobotac  is  a 
significant  aid  for  this,  especially  with  the  concentric 
annular  duct,  where  exit  flow  interaction  occurs  shortly 
past  the  duct  exit.  Once  this  transition  point  is  located, 
the  selector  valve  is  turned  to  allow  oil  to  flow  into  the 
weighing  tank.  The  mass  flow  rate  is  then  determined,  and 
Re.,  calculated  as. 


Re.-  =  mD/A  p 


/V 


The  hydrodynamic  entrance  length  was  determined 


graphically  from  the  pressure  vs.  position  plots  for  each 
duct.  In  fully  developed,  laminar  flow,  dP/dx  is  constant. 
In  the  hydrodynamic  entrance  region,  where  the  velocity 
profile  is  still  developing,  the  pressure  drop  is 
increasingly  greater  toward  the  entrance  and  thus, 
non-linear.  Where  the  non-linear  portion  of  the  pressure 
drop  curve  becomes  linear  defines  X  (see  Figure  11),  and 
knowing  the  duct  D*,  or  D^.j,  and  Re,  L*  can  be  determined. 


Figure  11.  Pressure  Drop  for  Fully  Developed,  Laminar  Flow 
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V.  Results  and  Discussion 

Friction  Factor 

Friction  factors  for  each  duct  can  be  determined  from 
the  pressure  drop  in  the  fully  developed,  laminar  flow 
region.  Friction  factor  data  was  used  to  verify  exact 
solutions,  based  on  Dh.  Friction  factors  were  determined 
from, 

4C,  =  (  -dP/dx )  (  2gD/P  V2) 

where  D  is  Dh  or  D*,i,  depending  on  the  application. 

Friction  factor  data  from  the  circular,  square,  and 
concentric  annular  ducts  are  compared  to  the  exact 
solutions,  based  on  Dh,  in  Figures  11,  12,  and  13, 
respectively.  Pressure  and  mass  flow  data  are  presented  in 
Appendix  C,  for  each  duct.  The  experimental  Cf's  agree 
very  well  with  the  exact  solutions  for  all  three  ducts. 

This  is  no  surprise,  since  experimental  C»  data  that 
verifies  the  exact  solutions  is  abundant  in  the  literature. 
More  importantly,  however,  this  data  provides  a  high  degree 
of  confidence  in  the  data  taking  and  reduction  techniques, 
which  carries  over  to  the  Re*,  and  L*  investigations. 

Since  the  C,  values,  based  on  Dh,  match  the  analytic  exact 
solutions,  by  definition,  C»,  based  on  Dh.i,  for  the  square 
and  concentric  annular  duct  will  equal  the  circular  duct 
value  (within  the  Dh.i  limitations  discussed  in  section 


CIRCULAR  DUCT 


Friction  Factors  for  Fully  Developed,  Laminar  Flow 


SQUARE  DUCT 


Friction  Factors  for  Fully  Developed,  Laminar 


CONCENT  QNNUL  DUCT 


Annular  Friction  Factors  for  Fully  Developed,  Laminar  Flow 


Transition  Reynolds  Number 

Transition  from  laminar  to  turbulent  flow  can  be 
precicted,  or  correlated,  by  a  transition  Reynolds  number. 
Re*, .  The  value  for  Re«,  for  circular  ducts  is 
approximately  2300.  This  is  only  under  normal  conditions, 
however.  If  disturbances  in  the  duct  entrance  region,  as 
well  as  the  approach  to  it,  are  minimized.  Re*,  can  be 
increased  to  significantly  higher  values,  the  upper  limit 
of  which  has  not  been  established.  V.W.  Ekman,  for 
example,  reached  a  Re»,  of  40,000  with  an  exceptionally 
disturbance  free  entrance  (15).  The  lower  bound  is 
approximately  2000,  below  which  laminar  flow  is  maintained 
even  with  very  strong  disturbances  present  (2:433) 

This  author  has  found  no  data  on  Re*,  for  non-circular 
ducts.  If  the  Re*,,  based  on  D.,,  is,  in  fact,  dissimilar 
for  different  duct  shapes,  a  desirable  trait  of  Dh.i  would 
be  to  correlate  the  Re.,  for  all  duct  shapes  to  the 
circular  duct  value  of  2000. 

The  three  ducts  each  have  abrupt  entrances  which 
induce  a  large  disturbance  at  the  entrance.  With  such  a 
disturbance,  transition  should  take  place  at  the  minimum 
Re.,  for  each  duct.  For  example,  historical  data  indicates 
that  this  will  occur  at  Re.,  =  2000  for  a  circular  duct. 

The  abrupt  entrance  should  eliminate  the  Re.,  dependence  on 
duct  roughness  . 


Transition  Reynolds  number,  based  on  Dh.i,  is 
defined  as: 

Re,„*.  =  (  P  VD„/  h  )  (Dh^/D*) 

Transition  Reynolds  number  were  first  calculated  based  on 
D„  to  compare  to  the  circular  duct  Re*,.  This  was  done 
because  this  author  found  no  data  on  Re*,  for  non-circular 
ducts.  The  Dh.i  correction  was  then  applied  to  try  to 
correlate  the  data  to  the  circular  duct  value.  For  these 
calculations,  ^  was  taken  from  Table  VI  for  all  three 
ducts.  This  was  done  after  the  extrapolated  values  from 
Table  VI  were  verified  from  the  circular  duct  data,  using 
the  Hagen-Poiseui lie  law  (see  Appendix  D) . 

Ten  to  15  data  runs  at  transition  were  made  for  each 
of  the  three  ducts,  with  the  average  Re*,,  based  on  both  D* 
and  D*,!,  shorn  in  Table  VII.  Individual  Re»,  values  can 
be  found  in  Appendix  C.  The  circular  duct  Re*,  was  2028, 
which  matches  the  historical  value  extremely  well.  The 
square  and  concentric  annular  ducts  had  Re*,s,  based  on  Dh, 
approximately  12.5%  below  the  circular  duct  Re*,.  What  is 
interesting  to  note  here  is  that  the  square  and  concentric 
annular  ducts,  each  have  similiar  Re*,,  even  though  they 
are  significantly  different  duct  shapes.  Recall,  these  two 
ducts  were  chosen  because  they  represented  two  extremes 
when  considering  C,Re,  based  on  Dh.  The  new  hydraulic 
diameter  appears  to  work  well  correlating  the  square  duct 


TABLE  VII 


Transition  Reynolds  Numbers 


Duct 

Re»r  -  Do. 

Retr  —  Dh.l 

Circle 

2028 

2028 

Square 

1750 

1925 

Concentr ic 
Annulus 

1733 

1165 

Re%r  to  that  of  the  circular  duct.  However,  for  the 
concentric  annular  duct,  the  Re„  diverges  from  the 
circular  duct  value.  This  data  indicates  that  Ret,,  based 


on  Dh, 


remains  the  best  choice. 


Hydrodynamic  entrance  length  is  defined  as  the  duct 
length  required  to  achieve  a  duct  centerline  velocity  99% 
of  the  fully  developed  value.  The  non-dimensional  form  is 

L*  =  ( X/D ) /Re 

where  D  is  either  Dh/  or  D*,,»,  depending  on  the 
application.  Hydrodynamic  entrance  length  is  important  in 
that  it  represents  a  position  where  fully  developed  flow 
can  be  assumed.  For  a  circular  duct,  Liu  (16)  calculated 
L*  =  .0541,  Heaton  (17)  predicted  L*  =  .0575,  and  Langhaar 
(18)  calculated  L*  =  .058.  For  equilateral  triangular 
ducts,  Shah  (19)  reported  L*  =  .0398.  Table  VIII  and  Table 
IX  present  analytic  L*  for  rectangular  and  concentric 
annular  ducts,  respectively.  Prom  these  data  it  is  easily 
seen  that  L*  varies  considerably  from  duct  shape  to  shape. 
In  fact,  as  Table  VIII  and  Table  IX  show,  authors  often 
disagree  on  L*  for  a  particular  duct  geometry.  Data  in 
Table  VIII  and  Table  IX  represent  analytic  results. 

Little  experimental  data  exists  to  verify  the  most  accurate 
analytic  method.  This  demonstrates  the  need  for  a 
parameter  such  as  D„.i  which  can  correlate  data  for 
different  duct  shapes  to  the  circular  duct  value. 

For  each  of  the  three  ducts,  10  to  15  laminar  flow 
data  runs  were  made,  at  various  Re,  to  determine  L*  (see 
Appendix  E ) . 


TABLE  VIII. 

Rectangular  Duct  L*  for  Fully  Developed  Laminar  Flow 


L* 


1/ a 

Wiginton 

Dalton 

Fleming 

Sparrow 

Han 

(21) 

(22) 

(23) 

1 

0.090 

0 . 0752 

0 .750 

- 

- 

0 . 0735 

0.500 

0.085 

0 . 095 

0 . 0660 

0.250 

0.075 

- 

0.0427 

0.200 

0.080 

0.080 

- 

0.125 

- 

- 

0.0227 

0 

— 

— 

0.0099 

I 


TABLE  IX. 

Concentric  Annular  Duct 
L*  for  Fully  Developed  Laminar  Flow 


ri/ri 

Liu 

(16) 

Manohar 

(25) 

u  _ 

Roy 

(26) 

Heaton 

(17) 

Coney 

(27) 

Sparro 

(28) 

0 

.0541 

- 

.0575 

- 

0 . 001 

- 

- 

- 

.0296 

- 

.0375 

0.01 

- 

- 

- 

- 

- 

.  0303 

0 .02 

- 

- 

- 

.0206 

- 

- 

0.05 

.  0206 

- 

- 

.0172 

.  0329 

.0241 

0.10 

.0175 

.0164 

.0180 

.0146 

.0253 

.0210 

0.20 

- 

- 

.0158 

- 

.0214 

.0171 

0.25 

- 

- 

- 

.0118 

.0204 

- 

0 . 30 

- 

.  0122 

.0140 

- 

.0194 

- 

0 . 40 

- 

- 

.0128 

- 

.  0178 

.0131 

0.50 

.  0116 

.0110 

.0121 

.  0103 

.  0168 

- 

0 . 60 

- 

- 

.0114 

- 

.  0161 

- 

0.70 

- 

.0103 

- 

- 

.  0156 

- 

0 .75 

.  0109 

- 

- 

- 

.0152 

- 

0 .80 

_ 

_ 

.0113 

.0150 

.0118 

1.00  .0108 


.  0099 


.0147 


TABLE  X 


Experimental 

Hydrodynamic 

Entrance  Length 

Duct 

L*  -  Dh. 

L*  -  D ^ 

Circle 

.0574 

.0574 

Square 

.0917 

.0758 

Concentr  ic 
Annulus 

( .011  -  Anal. 

)  ( .0163  -  Anal 

The  experimental  L*  values,  based  on  Dh  and  Dh.i,  are 
shown  in  Table  X.  The  hydrodynamic  entrance  length  for  the 
circular  and  square  ducts  were  easily  determined  from  plots 
like  those  shown  in  Figures  15  and  16,  respectively.  The 
concentric  annular  duct,  however,  posed  some  problems.  By 
examining  Figure  17,  it  is  seen  that  the  pressure  drop 
stays  constant,  for  the  concentric  annulur  duct,  from 
station  1  to  6,  even  at  high  Re.  A  hydrodynamic  developing 
region  is  indeterminable  from  the  data,  thus  making  it 
impossible  to  graphically  determine  L*.  By  examining  the 
analytic  L*  data  from  Table  IX,  a  possible  explanation  is 
found.  For  ri/r2  =  .84,  L*  is  approximately  .011.  For  a 
Re  of  1700,  just  before  transition,  X,  for  this  duct,  is 
approximately  4.6  inches.  Recall  the  first  two  pressure 
taps  were  located  at  3.0  and  5.3  inches  downstream  of  the 
entrance.  The  flow  is  fully  developed  before  it  reaches 
the  second  pressure  tap!  To  graphically  determine  L*,  with 
only  the  first  pressure  tap  within  the  hydrodynamic 
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Figure  16.  Hydrodynamic  Entrance  Length  for  Square  Duct 


Annular  Duct 


entrance  region,  and  close  to  the  end,  at  that,  is 
virtually  impossible.  The  first  three  pressure  taps  were 
located  as  close  to  the  entrance  as  physically  possible,  so 
a  L*  cannot  be  graphically  determined  for  this  particular 
duct . 

The  hydrodynamic  entrance  length,  based  on  Dh,  for  the 
circular  duct  was  0.0574  which  indicates  that  Heaton  (17) 
and  Langhaar  (18)  have  good  analytic  techniques  in 
predicting  L* .  The  square  duct  had  a  L*  of  0.0917.  This 
value  indicates  that  Wiginton  and  Dalton  (21)  have  the  best 
analytic  technique  for  determining  L*  for  square  ducts. 

Using  Dh.i  as  the  characteristic  dimension  instead  of 
Dh,  does  not  correlate  the  data  to  the  circular  duct  value. 
The  correlation  trend,  however,  is  favorable.  The  square 
duct  L*  is  reduced  toward  the  circular  duct  value,  while 
the  concentric  annular  duct  L*  is  increased. 

Since  the  correlation  of  L*,  based  on  Dh.i,  produced 
unsatisfactory  results,  perhaps  a  different  approach  to  the 
L*  problem  is  needed.  Consider  flow  through  a  circular 
duct  as  shown  in  Figure  18.  Here  the  boundary  layer 
growth,  6  ,  is  generally  described  as, 

6  =  C(i^t)-S  (2:83)  (10) 

where  C  is  a  function  of  the  pressure  gradient,  dP/dx,  and 
the  transverse  velocity,  v.  We  will  let  time,  t,  be  equal 
to  x/V,  for  this  problem.  (It  should  be  noted  that  V  is 
not  the  actual  velocity  that  the  boundary  sees  from 
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Figure  18.  Boundary  Layer  Growth  in  the  Entrance  Region 


the  entrance  to  position  X.  The  velocity  actually  varies. 


non-linearly,  from  V,  at  the  entrance,  to  2V,  at  the 


centerline,  for  a  circular  duct,  but  we  are  using  V  here 


for  convenience  and  to  obtain  "ballpark"  expressions  for 


L* . )  The  hydrodynamic  entrance  length,  X,  becomes 


approximately  the  position  where  the  boundary  layer  meets 


the  duct  centerline,  <5  =  D/2.  Equation  (10)  becomes. 


D/2  =  CUX/V)-9 


which  can  be  rewritten. 


X/D  =  ( VD/ v ) ( 1/ 4C? ) 


(11) 


where  X  OC  DJ . 


L*  =  ( X/D ) /Re  =  1/4C1 


(  12  ) 


For  the  circular  duct,  L*  =  .0574,  and  C  is  approximately 


2.1. 


Now  let  us  keep  C  constant  at  2.1,  and  examine  L*  for 


non-circular  ducts,  where  6  becomes  the  average  distance  a 


disturbance  must  travel  to  reach  the  duct  centerline.  Let 
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Figure  19.  A  Square  Duct 

us  examine  a  square  duct,  shown  in  Figure  19,  for  instance, 
where,  from  the  corners,  a  disturbance  must  travel 
(2-9)(Dh/2)  to  reach  the  centerline.  For  non-circular 
ducts,  X  is  proportional  to  D*2,  so, 

L*  =  (2) ( .0574)  =  .1148 

This  would  be  the  predicted  L*  if  all  the  disturbances  in  a 
square  duct,  travelled  ( 2-")  (D*/2 )  )  to  the  center. 

Obviously  this  is  not  the  case.  Those  in  the  corners 
travel  (2-*)(Dh/2),  but  those  on  the  sides  only  travel 
(Dn/2)  to  the  centerline.  So,  perhaps,  the  average 
distance  a  disturbance  travels  in  a  square  duct  is 
somewhere  between  these  two  values,  say  at  ( 1 . 5  )  •*(  D„/2  ) , 
for  instance.  If  this  is  the  case, 

L*  =  (1.5) ( .0574)  =  .0861 

This  L*  is  very  close  to  our  experimental  value,  L*  = 

.0917!  Likewise,  for  an  infinite  parallel  plate  duct,  the 
average  disturbance  travel  distance  is  .  5(D„/2),  so. 


L 


.  01435 


This  also  is  quite  close  to  the  analytically  predicted  L* 


So,  perhaps,  a  better  way  to  correlate  L*  data  would 
be  to  scale  the  constant,  C,  with  a  new  hydraulic  diameter 
that  estimates  the  average  boundary  layer  growth  to  the 
centerline  of  non-circular  ducts. 
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VI  . 


Conclusions  and  Recommendations 


Conclusions 

In  this  investigation,  a  new  hydraulic  diameter,  Dh.i, 
that  better  represented  the  mean  center  to  surface  distance 
for  laminar  flow  in  non-circular  ducts  was  empirically 
determined.  Experimentation  was  performed  to  examine  the 
effect  on  Re«r  and  L*  of  replacing  the  traditional 
hydraulic  diameter,  D*,  with  Dh,>.  The  experiment  was 
performed  using  a  circular  duct,  square  duct,  and 
concentric  annular  duct  that  simulated  flow  between 
infinite  parallel  plates. 

The  new  hydraulic  diameter  equation  correlated  C,Re 
data  to  the  circular  duct  value  for  many  non-circular  duct 
shapes.  All  the  shapes  examined  correlate  within  2.4 
percent  of  the  circular  duct  value,  with  most  within  1.5 
percent . 

Experimental  Re*,  were  obtained  for  the  circular,  square, 
and  concentric  annular  ducts.  The  circular  duct  Re.,, 
based  on  Dh,  was  very  close  to  2000  with  the  square  and 
concentric  annulus  Re*,  approximately  12.5  percent  lower 
than  this  value.  Square  duct  Re*,,  based  on  Dh.i, 
correlated  very  close  to  the  circular  duct,  but  the 
concentric  annulus  Re*,  did  not.  Concentric  annulus  Re*, 
actually  diverged  from  the  circular  duct  Re*,.  Based  on 
this  data,  it  appears  that  Dhli,  which  correlates  C,Re  data 


very  well,  is  not  a  good  character  is  t  ic  dimension  for  Re*,. 
Perhaps  an  entirely  different  hydraulic  diameter  is  needed 
if  Re*r  is  to  be  correlated  to  the  circular  duct  value. 

Hydrodynamic  entrance  lengths  were  determined  for  only 
the  circular  and  square  ducts.  The  concentric  annular  duct 
flow  became  fully  developed  before  the  second  static 
pressure  tap,  making  accurate  prediction  of  L*  almost 
impossible.  The  analytic  L4  value  for  the  concentric 
annulus  was  thus  used  for  D»,,i  correlations.  Hydrodynamic 
entrance  lenghts,  based  on  Dt.P»,  did  not  correlate  well  to 
the  circular  duct  value.  The  correlation  did,  however, 
display  the  correct  trend  toward  the  circular  duct  value, 
and,  perhaps,  multiplication  of  D*,i  by  some  constant  might 
bring  the  L*  values  in  line. 

The  new  hydraulic  diameter  is  quite  accurate  in 
correlating  C,Re  data,  but  cannot  be  relied  upon  to 
correlate  Re*,  and  L*  data.  Transition  Reynolds  numbers, 
based  on  Dh,  were  determined  for  the  ducts,  and  will  be 
useful,  since  no  data  for  non-circular  duct.  Re*,,  has  been 
reported.  Likewise,  experimental  circular  and  square  duct 
L*  values,  based  on  D*,  will  add  to  the  limited  data  base 
available . 

Recommendations 

The  new  hydraulic  diameter  did  not  correlate  Re*,  and 
L4  data.  Separate  hydraulic  diameters  are  needed  to 
correlate  both  Re*,  and  L4  data.  The  L4  correlation 


appears  to  require  only  a  multiplier  applied  to  the  Dh,i 
definition. 

Instead  of  using  Dh,i  directly  in  the  L*  correlation, 
perhaps  determining  a  D*,i  that  correlated  the  constant,  C, 
in  the  boundary  layer  growth  equation  (Equation  (10)), would 
yeild  better  results,  and  should  be  investigated  further. 

Pressure  gages,  located  close  to  the  static  pressure 
taps,  should  be  installed,  instead  of  relying  on  the 
10-foot  pressure  lines  that  currently  lead  to  the  manometer 
board.  This  would  provide  better  response  indications  to 
the  pressure  fluctuations  which  indicate  the  onset  of 
transition.  In  addition,  a  static  pressure  tap  should  be 
added  to  each  duct  at  the  exit.  This  would  aid  in  the 
detection  of  flow  transition. 

Static  pressure  taps  should  be  installed  very  close  to 
the  entrance  of  the  concentric  annular  duct  to  help  in  the 
L*  determination.  This  will  require  a  major  material 
overhaul . 

Unusual  duct  shapes,  such  as  triangular,  hexagonal,  or 
semi-circular,  should  be  used  in  the  investigation  of  CfRe, 
Re»,,  and  L*,  to  increase  the  experimental  data  base  used 
in  the  D*.i  correlations.  It  would  be  interesting  to  3ee 
where  the  data  for  these  ducts  fall,  compared  to  that  of 
the  square  and  concentric  annulus. 
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Appendix  A. 

New  Hydraulic  Diameter  Calculations 


Notation : 

A  =  Duct  area 

P  =  Duct  perimeter 

Ac  =  Area  of  the  modelling  circle 

Pc  =  Perimeter  of  the  modelling  circle 

Ae  =  A  -  Ac 

Pe  =  P  +  Pc 

W  =  Width  of  infinite  parallel  plates  (approaches  infinity) 

TABLE  XI . 

Compact  Duct  Hydraulic  Diameters 
Shape  A  P  Ac  Pc  Ae  Pe  Ae/Pe  Dk..  D* 

Circle  .78S2  3.14S  .78S2  3.14S  .00S2  6.28S  .000S  1.00S  1.00S 

Square  1.0S2  4.00S  .78S2  3.14S  .21S2  7.14S  .031S  1.10S  1.00S 

Triang  .43S2  3.00S  . 26S2  1.81S  .17S2  4.81S  .036S  .693S  .577S 

Hexag  2.6S2  6.00S  2.4S2  5.44S  .24S2  11. 4S  .021S  1.80S  1.73S 


Semi-  . 40S2  2.57S  .26S2  1.80S  .14S2  4.37S  .033S  .633S  .611S 
Circle 


Infinite  A  =  WS  P  =  2W  +  2S 

Parallel 

Plates  Ac  =  .79WS  Pc  =  2W 


Ae  =  . 2 1WS  Pe  =  4W  +  2S  Ae/Pe  =  . 21WS/ ( 4W+2S ) 
As  W  goes  to  infinity,  Ae/Pe  reduces  to  .05365S 
Dh,.  =  1.34S  D„  =  2S 
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TABLE  XII 


Rectangular  Hydraulic  Diameters 


g_  a  E. 


1.00 

l.OOS* 

4. OS 

1.25 

1.25S* 

4.5S 

1.50 

1.50S* 

5. OS 

1.75 

1.75S* 

5.5S 

2.00 

2. OOS* 

6.  OS 

2.25 

2.25S* 

6.5S 

2.50 

2.50S* 

7.  OS 

2.75 

2.75S* 

7.5S 

3.00 

3. OOS* 

8  .OS 

3.25 

3.25S* 

8.5S 

3.50 

3.50S* 

9.  OS 

3.75 

3.75S* 

9 . 5S 

4.00 

4  .  OOS  * 

10. OS 

4.25 

4.25S* 

10. 5S 

4.50 

4. SOS* 

11. OS 

4.75 

4.75S* 

11. 5S 

5.00 

5. OOS  * 

12.  OS 

5.25 

5.25S  * 

12. 5S 

5.50 

5.50S  * 

13. OS 

5.75 

5.75S  * 

13. 5S 

6.00 

6. OOS  * 

14. OS 

6.25 

6.25S  * 

14. 5S 

6.50 

6.50S  * 

15. OS 

6.75 

6.75S  • 

15. 5S 

7.00 

7. OOS  * 

16.  OS 

7.25 

7.25S  * 

16. 5S 

7.50 

7.50S  * 

17. OS 

7.75 

7.75S  * 

17. 5S 

8.00 

8. OOS* 

18.  OS 

8.25 

8.25S  * 

18. 5S 

8.50 

8.50S* 

19.  OS 

8.75 

8.75S* 

19. 5S 

9.00 

9. OOS* 

20. OS 

9.25 

9.25S* 

20. 5S 

9.50 

9.50S* 

21. OS 

9.75 

9 . 75S  * 

21. 5S 

10.00 

10. OOS* 

22. OS 

25.00 

25. OOS* 

52. OS 

50.00 

50. OOS* 

102 .OS 

100.0 

100. OS* 

202. OS 

Ac 

Pc 

&e 

0.79S* 

3  .OS 

0. 21S* 

0.98S* 

3.4S 

0.27S* 

1.18S* 

3.9S 

0.32S* 

1.37S* 

4. 3S 

0.38S* 

1.57S* 

4.8S 

0. 43S* 

1.77S* 

5.3S 

0. 48S* 

1.96S* 

5.7S 

0. 54S* 

2.16S* 

6.2S 

0.59S* 

2.36S* 

6 . 7S 

0.64S* 

2.55S* 

7. 2S 

0 . 70S  * 

2.75S* 

7.7S 

0.75S* 

2.95S* 

8.2S 

0 . 80S  * 

3.14S* 

8.7S 

0.86S* 

3.34S* 

9.2S 

0.91S* 

3.53S* 

9.7S 

0.97S* 

3.73S* 

10. 2S 

1.02S* 

3.93S* 

10. 7S 

1.07S* 

4 . 12S* 

11. 2S 

1.13S* 

4 . 32S* 

11. 6S 

1.18S* 

4.52S* 

12. IS 

1.23S* 

4.71S* 

12. 6S 

1.29S* 

4.91S* 

13. IS 

1.34S* 

5. US* 

13. 6S 

1.39S* 

5.30S* 

14. OS 

1. 45S* 

5.50S* 

14. 5S 

1.50S* 

5.69S* 

15. OS 

1.56S* 

5.89S* 

15. 5S 

1.61S* 

6.09S* 

15. 9S 

1.66S* 

6.28S* 

16. 4S 

1.72S* 

6.48S* 

16. 8S 

1.77S* 

6.68S* 

17. 3S 

1. 82S* 

6.87S* 

17. 8S 

1. 88S* 

7.07S* 

18. 2S 

1.93S* 

7.26S* 

18. 7S 

1.99S* 

7 . 46S  * 

19. IS 

2.04S* 

7.66S* 

19. 6S 

2.09S* 

7.85S* 

20. OS 

2.15S* 

9 . 64S* 

50. OS 

5. 37S* 

39.3S* 

100  .  OS 

10. 7S* 

78. 5S* 

200. OS 

21.  SS* 

Pe  Ae/Pe  D».i 

7. OS  . 031S  1.10S 

7.95  .  034S  1.24S 

8 . 95  . 036S  1.29S 

9.8S  . 038S  1.33S 

10. 85  . 040S  1.35S 

11. 85  . 041S  1.37S 

12. 75  . 042S  1.38S 

13. 75  . 043S  1.39S 

14. 75  . 044S  1.39S 

15. 75  . 044S  1.40S 

16. 75  . 045S  1.40S 

17. 75  . 045S  1.40S 

18. 75  . 046S  1.40S 

19. 75  . 046S  1.40S 

20. 75  . 047S  1.40S 

21. 75  . 047S  1.41S 

22. 75  . 047S  1.41S 

23. 75  . 048S  1.40S 

24. 65  . 048S  1.40S 

25. 65  . 048S  1.40S 

26. 65  . 048S  1.40S 

27. 65  . 049S  1.40S 

28. 65  . 049S  1.40S 

29. 55  . 049S  1.40S 

30. 55  . 049S  1.40S 

31. 55  . 049S  1.40S 

32. 55  . 050S  1.40S 

33. 45  . 050S  1.40S 

34. 45  . 050S  1.40S 

35. 35  . 050S  1.40S 

36. 35  . 050S  1.40S 

37. 35  . 050S  1.40S 

38. 25  . 051S  1.40S 

39. 25  . 051S  1.40S 

40. 15  . 051S  1.40S 

41. 15  . 051S  1.40S 
42. OS  . 051S  1.40S 

102. OS  . 053S  1.37S 
202. OS  . 053S  1.36S 
402. OS  . 053S  1.35S 


D* 

1.00S 
1.11S 
1.20S 
1.27S 
1.33S 
1.38S 
1.43S 
1.47S 
1.50S 
1.53S 
1.56S 
1.58S 
1 . 60S 
1.62S 
1.64S 
1.65S 
1.67S 
1.68S 
1.69S 
1 . 70S 
1.71S 
1.72S 
1.73S 
1.74S 
1.75S 
1.76S 
1.76S 
1.77S 
1.78S 
1.78S 
1.79S 
1.79S 
1 . 80S 
1 . 80S 
1 .81S 
1.81S 
1 .82S 
1.92S 
1.96S 
1 .98S 


Specific  Gravity  Measurements 


TABLE  XIII. 

AFIT  SAE  10W10  Oil  Specific  Gravity  Measurements 
Measurement  Water  (grams)  Oil  (grams)  Sp . 


Data  Reduction 


TABLE  XIV. 

Data  for  Circular  Duct,  Laminar  Flow 
Run  Number 

1 _ 2 _ 3 _ 4 _ 5 _ 6 


Pi  -  in.  Hg 

3.20 

2.77 

2.29 

1.92 

2.14 

2.24 

P 2 

2.83 

2.49 

2.06 

1.78 

1.92 

2.00 

P  3 

2.53 

2.20 

1.83 

1.60 

1.72 

1.78 

P« 

2.25 

1.95 

1.64 

1.45 

1.50 

1.55 

P« 

1.97 

1.71 

1.45 

1.25 

1.35 

1.39 

P. 

1.78 

1.58 

1.35 

1.19 

1.28 

1.29 

Pt 

1.48 

1.32 

1.16 

1.01 

1.09 

1.10 

P. 

1.28 

1.16 

1.02 

0.97 

0.98 

0.99 

P, 

1.02 

0.92 

0.83 

0.78 

0.80 

0.80 

P.o 

0 . 76 

0.74 

0.70 

0.69 

0 . 70 

0 . 70 

Temp  -  F 

81.0 

83.0 

85.0 

87.5 

92.0 

93.0 

m  -  lbm/sec 

0.84 

0.74 

0.61 

0 . 49 

0 .61 

0.65 

V  -  ft/sec 

4.33 

3.82 

3 . 16 

2.53 

3.15 

3.35 

C, 

.0366 

.0380 

.  0440 

.0560 

.0405 

.0365 

RE 

1770 

1636 

1426 

1132 

1574 

1698 

TABLE  XV. 


Data  for  Circular  Duct,  Transition  Flow 
Data  Set  1 

Run  Number 


1 _ 2 


Pi  -  in.  Hg 

3.20 

2.89 

P» 

2.90 

2.60 

P* 

2.51 

2.27 

P  4 

2.22 

2.02 

Ps 

1.91 

1.75 

p* 

1.70 

1.60 

p? 

1.42 

1.32 

p. 

1.28 

1.19 

p* 

1.00 

0.91 

P.o 

0.80 

0.79 

Temp  -  F 

84.8 

86.7 

• 

m  -  lbm/sec 

0.88 

0.82 

V  -  ft/sec 

4 . 51 

4.23 

C, 

.0316 

.0330 

RE 

2023 

1958 

3 

4 

5 

6 

2.90 

2.80 

2.51 

2.30 

2 . 61 

2 . 50 

2.23 

2.04 

2.29 

2.22 

2.00 

1.81 

2.02 

1.96 

1.78 

1.60 

1.73 

1.73 

1.56 

1 . 40 

1.61 

1.55 

1.40 

1.30 

1.34 

1.30 

1.10 

1.11 

1.18 

1.15 

1.03 

1.01 

0.93 

0.92 

0.86 

0.83 

0.79 

0.74 

0.72 

0.71 

88.0 

90.0 

92.1 

94.8 

0.81 

0.80 

0.79 

0.74 

4.15 

4.13 

4.05 

3 .83 

0339 

.0324 

.0307 

.0295 

1962 

2014 

2027 

1976 

Fjrr\*  j  \jh  na  nji !  l*  r\*  rarnjirw  njnvTyr^ji'nai  rurrurr\»  >'\?r  Terror 


TABLE  XVIII 


Data  for  Square  Duct,  Laminar  Flow 
Data  Set  2 

Run  Humber 


1 

2 

3 

Pi  -in.  Hg 

1.71 

1.77 

3.50 

P* 

1.55 

1.60 

3.12 

Ps 

1.38 

1.41 

2.69 

P« 

1.31 

1.30 

2.39 

Ps 

1.13 

1.13 

2.04 

P. 

1.09 

1.08 

1.81 

Pr 

0.92 

0.90 

1.50 

p. 

0.86 

0.87 

1.30 

p. 

0.74 

0.73 

1.04 

Pio 

0.62 

0.62 

0.72 

Temp  -  F 

100  . 

100. 

78.0 

m  -  lbm/sec 

0.56 

0.57 

0.89 

V  -  ft/sec 

2.66 

2.71 

4 .23 

Cf 

.0373 

.0360 

.0389 

RE 

1522 

1549 

1493 
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Data  Set  1 


TABLE  XX. 


Data  for  Square  Duct,  Transition  Flow 
Data  Set  2 


Run  Number 


TABLE  XXI 


Data  for  Concentric  Annular  Duct,  Laminar  Flow 

Data  Set  1 

Run  Number 


1 

2 

3 

4 

5 

6 

in .  Hg 

28.58 

30.24 

19.48 

23.18 

20.94 

24.70 

27.80 

29 . 41 

19.00 

22.60 

20.44 

24.09 

26.61 

28.18 

18.32 

21.80 

19.75 

23.40 

19.23 

20.26 

13.11 

15.50 

14.15 

16.60 

11.12 

11.50 

8.30 

8 . 30 

7.99 

9.27 

3.79 

3.97 

2.23 

3.29 

2.90 

3.36 

-  F 

90.0 

92.0 

91.8 

93.5 

96.3 

88.0 

lbm/sec 

2.12 

2.33 

1.54 

1.86 

1.89 

1.80 

ft/sec 

10.2 

11.2 

7.42 

8.94 

9.10 

8.67 

.0639 

.0558 

.0853 

.0660 

.0594 

.0760 

L0. 


TABLE  XXII. 


Data  for  Concentric  Annular  Duct,  Laminar  Flow 

Data  Set  2 

Run  Number 


1 

2 

3 

4 

5 

6 

Pi  -  in.  Hg 

20.64 

22.39 

25.13 

28.38 

30.41 

31.40 

P» 

20.16 

21.85 

24.53 

27.65 

29.65 

30.65 

Pa 

19.50 

21.16 

23.74 

26.79 

28.65 

29.68 

P  * 

13.88 

15.10 

16.83 

19.00 

20.30 

21.08 

Pa 

8.18 

8.48 

9.44 

10.64 

11.75 

12.30 

P. 

2.87 

3.06 

3.40 

3.80 

4.04 

4.18 

Temp  -  F 

88.0 

88.2 

88.7 

89.2 

89.8 

88.0 

« 

m  -  lbm/sec 

1.56 

1.70 

1.92 

2.14 

2.31 

2.28 

V  -  ft/sec 

7.52 

8.19 

9.25 

10.3 

11.1 

11.0 

« 

*1 

.0842 

.0775 

.0677 

.0616 

.0569 

.0605 

RE 


1071  1170 


1332  1500  1627  1565 


TABLE  XXIII. 


.V. 


Data  for  Concentric  Annular  Duct,  Transition  Flow 

Data  Set  1 


Run  Number 


1 

2 

3 

4 

5 

P.  -  in.  Hg 

34.10 

33.30 

31.35 

30.02 

28.85 

P* 

33.16 

32.40 

30.45 

29.20 

28.10 

P* 

31.96 

31.25 

29.38 

28.18 

27.05 

P4 

23.15 

22.50 

21.24 

20.30 

19.50 

Ps 

13.08 

12.56 

12.22 

11.70 

11.28 

P. 

4.49 

4.38 

4.15 

3.99 

3.85 

Temp  -  F 

90.0 

90.2 

92.0 

93.0 

94.0 

m  -  lbm/sec 

2.54 

2.49 

2.42 

2.37 

2.33 

V  -  ft/sec 

12.2 

12.0 

11.7 

10.8 

11.2 

C, 

.0538 

.0543 

.0544 

.0606 

.0536 

RE 

1801 

1770 

1758 

1645 

1736 

TABLE  XXIV 
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Data  for  Concentric  Annular  Duct,  Transition  Flow 

Data  Set  2 

Run  Number 


1 

2 

3 

4 

5 

Pi  -  in.  Hg 

27.95 

24.96 

34.00 

33.20 

35.00 

P* 

27.20 

24.30 

33.06 

32.30 

34.05 

Pj 

26.28 

23.51 

31.96 

31.30 

33.00 

P« 

18.20 

16.66 

22.73 

22.20 

23.42 

Pa 

10.70 

9.46 

13.05 

12.90 

13.43 

P» 

3.70 

3.35 

4.46 

4.40 

4.60 

Temp  -  F 

95.0 

96.5 

90.0 

88.3 

89.5 

m  -  lbm/sec 

2.20 

2.08 

2.55 

2.43 

2.60 

V  -  ft/sec 

10.6 

10.0 

12.3 

11.7 

12.5 

C, 

.0557 

.0572 

.0523 

.0565 

.0518 

RE 

1660 

1649 

1808 

1676 

1828 
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Viscosity  Comparison 


This  appendix  compares  viscosity  calculated  using  the 
Hagen-Poiseuille  Law  verses  extrapolated  values  from  Table 
VI.  The  Hagen-Poiseuille  Law  is  used  with  circular  duct 
flow  data  at  various  temperatures.  Table  VI  values  are 
extrapolated  from  two  data  points  using  a 
Viscosity-Temperature  chart. 


Re 


1426 

.0770 

3 . 162 

2 .470 

2 . 510 

1574 

.0700 

3.146 

2 . 254 

2 . 264 

1636 

.  1000 

3 .817 

2 .654 

2.642 

1698 

.0700 

3.352 

2.120 

2 .240 

Since  the  extrapolated  values  from  Table  VI  predict 
those  calculated,  using  the  Hagen-Poiseuille  Law,  they  will 
be  used  in  all  other  calculations,  for  the  sake  of 


TABLE  XXV. 

Viscosity  Comparison 

dP/dx  V  n -  H-P  calcs  M -  Table  VI 

in . Ha/ft  ft/sec  lbf-sec/f t*  lbf -sec/ft1 

.0625  2.527  2.500x10-  2.420x10- 
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Hydrodynamic  Entrance  Length  Plots 
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l.V*  FOR  C'RCUl.ftR 


?0.  Hydrodynamic  Hn  trance  Length  for  Circular  Duct,  Re  =  1770. 
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Hydrodynami c  Entrance  Length  for  Circular  Duct,  Re  =  14?6. 
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re  23-  Hydrodynamic  Entrance  Length  for  Circular  Duct,  Re  =  I698. 
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Figure  29.  Hydrodynamic  Entrance  Length  for  Circular  Duct,  Re 
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Figure  30.  Hydrodynamic  Entrance  Length  for  Square  Duct,  Re 
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